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The thermal fluctuations that exist at very low temperature in disordered systems are often at- 
tributed to the existence of some two-level excitations. In this paper, we revisit this question via the 
explicit studies of the following ID models (i) a particle in ID random potentials (ii) the random field 
Ising chain with continuous disorder distribution. In both cases, we define precisely the 'two-level' 
excitations and their statistical properties, and we show that their contributions to various observ- 
ables are in full agreement at low temperature with the the rigorous results obtained independently. 
The statistical properties of these two-level excitations moreover yield simple identities at order T in 
temperature for some generating functions of thermal cumulants. For the random-field Ising chain, 
in the regime where the Imry-Ma length is large, we obtain that the specific heat is dominated 
by small non-universal excitations, that depend on the details of the disorder distribution, whereas 
the magnetic susceptibility and the Edwards- Anderson order parameter are dominated by universal 
large excitations, whose statistical properties only depend on the variance of the initial disorder via 
the Imry-Ma length. 

I. INTRODUCTION 
A. Disordered systems at very low temperature 

One of the most important feature of disordered systems is that they may present states that have an energy very 
close to the ground state energy but which are very different from the ground state in configuration space. For the 
spin-glasses, the debate between the droplet and replica theories concerns the probabilities and the properties of 
these states. In the droplet theory 0, the low-temperature physics is described in terms of rare regions with nearly 
degenerate excitations which appear with a probability that decays with a power-law of their size. In the replica 
theory 0, the replica symmetry breaking is interpreted as the presence of many pure states in the thermodynamic 
limit, i.e. the nearly degenerate ground states appear with a finite probability for arbitrary large size. 

More generally, the statistical properties of the nearly degenerate excitations (their numbers, their sizes, their 
geometric properties, the barriers separating them, etc..) are interesting in any disordered system, since they govern 
all properties at very low temperature. In particular, a linear behavior in temperature of the specific heat 

C{T)=bT + 0{T^) (1) 

seems rather generic for a large class of disordered models, ranging from spin-glasses where this behavior is measured 
experimentally Q and numerically to one-dimensional classical spin models where this behavior can be exactly 
computed via the Dyson-Schmidt method . This linear behavior of the specific heat has been also obtained recently 
for classical disordered elastic systems via the replica variational method ^] . In this respectjthe infinite-ranged spin 
glass models are anomalous since the first term of the specific heat is only of order 0, Q . The leading term of 
the specific heat at low temperature has also been studied recently in various quantum disordered models, and is still 
under debate between linear C ~ T .gj , quadratic C 0| and cubic C ~ [TH . 

The linear behavior of the specific heat is of course reminiscent of the same behavior in ordinary glasses, where 
the usual description involves a phenomenological theory of two- level systems 0, 0, ^3 ■ Whereas the identification 
of these two-level systems in a microscopic theory has remained problematic for structural glasses (see 0| for recent 
discussions), this problem seems simpler in disordered models, where concrete proposals have been made. In the 
context of spin-glasses for instance , it has been proposed a long time ago that the two-level excitations are clusters 
of spins that can be flipped with respect to the ground state, and that the associated joint probability iV(e, v) of the 
energy e and barrier v determines various static and dynamical properties at very low temperature 15]. However, the 
computation of this probability distribution N{e,v) is difficult except for very small clusters containing only a few 
spins • This is why the statistics of large excitations relies on some scaling assumptions in the droplet theory . 
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In recent years, there has been a lot of numerical efforts to characterize the distribution and the topology of these 
large low-energy excitations p^] . 

There exists a special class of disordered models for which exact remarkable identities for thermal fluctuations have 
been obtained The Hamiltonian of these models have a deterministic part which consists in quadratic interactions 
and a random part whose statistics is translation invariant. The simplest of these systems is the so called toy model 
[Tsj l defined by the one-dimensional Hamiltonian 

Htoy = fa;' + V{x) (2) 

where the random potential V{x) is a Brownian motion presenting the correlations 



{V{x)-V{v)f = 2a\x-v\ (3) 

The statistical tilt symmetry leads to the following result for the generating function of thermal cumulants averaged 
over the disorder [l3 
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In < e-^^ > = T— (4) 

This identity is particularly interesting at very low temperature, since it predicts a linear behavior in T of the second 
cumulant. In a previous work we have explicitly shown that the rare configurations with two nearly degenerate 
minima Ai? ~ T actually give the full exact second cumulant of I0J. The statistical tilt symmetry is also present in 
the models of directed polymers in random media, and we refer the reader to the references [23, |2l| for a detailed 
discussion of the identities on thermal fluctuations and their interpretation in terms of nearly degenerate paths. 



B. Goal and results 



The above presentation on the properties of disordered systems at very low temperature, although very incomplete, 
already shows that generically, various observables exhibit variations at order T in temperature, that can be related to 
the statistical properties of nearly degenerate excitations. The aim of this paper is to revisit this question in some one- 
dimensional models to identify precisely the 'two-level' excitations and their statistical properties, in order to compare 
their contribution to various observables at order T with the rigorous results that can be obtained independently. 



1. Results for one particle in a random potential 



We will first consider the equilibrium of a particle in random potentials of the following form 

H{x) = Ho{x) + V{x) (5) 

where i?o is the non-random part of the Hamiltonian, and where the random potential V{x) is a Brownian potential 
0. As explained above, the toy model i/p — f is very special since it satisfies exact identities Q as a consequence 
of the statistical tilt symmetry. Here we show that for a large class of deterministic part Ho, the following temperature 
expansion holds for the generating function of thermal cumulants 

\2 r + co 

ln<e-^-> =-Xx— + T— y''D{y) + 0{T^) (6) 

^ Jo 

where Xmin is the disorder average of the position Xmin where H{x) is minimum, and where D{y) is the probability 
to have two degenerate minima in the system separated by a distance y. The thermal fluctuations of these models are 
thus directly related to the presence of metastable states in rare disordered samples. In particular, the susceptibility 

X=^{<x' > - <x>^) (7) 
has for average a finite value at zero-temperature 



X= / y^D{y) + 0{T) 
Jo 



(8) 
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but only the samples with two nearly degenerate minima contribute to this zero-temperature value. The typical 
samples with only one minimum have a susceptibility that vanishes at zero temperature. This type of behavior was 
found numerically for the random directed polymer model [23| • We will show that the relation 10 is in full agreement 
with the exact results corresponding to the two soluble cases 

(i) the pure Brownian potential Hq = on a finite interval x € [0, L]; 

(ii) the biased Brownian potential Hq = fx on the semi- infinite line x G [0, +oo[. 



2. Results for the random field Ising chain 



We will then consider the one-dimensional random field Ising model, in the regime where the Imry-Ma length Lim 
representing the typical domain size at zero temperature is large. In contrast with models ^ with only one degree of 
freedom, where the low-temperature properties are governed by the rare samples having two nearly degenerate ground 
states, the systems with many degrees of freedom are well described by disorder averaged values in the thermodynamic 
limit, since the average over the disorder is like a spatial average. As a consequence, the low temperature properties 
are governed by the rare regions that presents nearly degenerate excitations. We describe in detail the rare nearly 
degenerate excitations and evaluate their contribution to the low temperature behavior of various observables. The 
density p{E = 0, 1) of excitations of length I contains two contributions, involving either a single domain wall, that 
has two nearly degenerate optimal positions, or a pair of neighboring domain walls, than can appear or disappear 
with almost no energy cost. Our analysis allows to obtain that the specific heat |0] 

C{T) = T^J dlp{E = 0, + 0{T^) (9) 

is dominated by small non-universal excitations of length 1^1, that depend on the details of the disorder distribution, 
whereas the Edwards- Anderson order parameter 

r+oo 

qEA = <Si>^ = l-2T dllp{E = 0,l) + O{T'^) (10) 

^0 

and the magnetic susceptibility 

X= — (^< m2 > - < m >2j = 2y dl Pp{E ^0,l) + O{T) (11) 

are dominated by large excitations whose length I is of order of the Imry-Ma length Ljm, and whose properties are 
universal with respect to the initial disorder distribution, since they only depend upon its variance. We will then 
compare with the exact results available from the Dyson-Schmidt method . 



3. Organization of the paper 

The paper is organized as follows. The Section ^ is devoted to the equilibrium at low temperature in the one- 
dimensional potentials ((SJ). The Section imi concerns the analysis of two-level excitations in the random field Ising 
chain, and their influence on various observables. The Appendix ^ contains some exact results on thermal cumulants 
in the Brownian potential, that are useful to show the exactness of the two- level description presented in Section HII 



II. EQUILIBRIUM OF A PARTICLE IN ONE-DIMENSIONAL RANDOM POTENTIALS 

A. Effective model at order T in temperature 

In typical samples with Hamiltonian ||SJ| , there is only one minimum. In , it was moreover shown that the typical 
extension of the Boltzmann distribution around a minimum for a pure Brownian potential is given by the thermal 
length It = that comes from the Boltzmann factor e^^^'^-* with the typical behavior V{x) ^ \fax. For the 
other models with regular deterministic part, the leading behavior at small distance is still given by the Brownian 
potential y(x\ and thus the same argument can be applied. As a consequence, for all potentials (0), the typical 
fluctuations around a minimum are of order Aa; ~ at low temperature, and thus do not contribute at order T . 
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FIG. 1: Example of a random potential presenting two nearly degenerate minima situated at the positions xi and X2- The 
low temperature equilibrium properties are completely determined by the probability D(y) over the samples that there exist 
two nearly degenerate minima separated by a distance y. 



We now consider the configurations presenting two almost degenerate minima H{xi) ~ H{x2) situated at Xi and 
X2 '■ see Figure ^ To compute their contribution at order T, we may neglect the spatial extension of order It ~ T'^ 
around each minimum, and replace the thermal measure by two delta peaks 

e-/3'»e// ^ pS(^x - xi) + (1 - p)6{x ~ X2) (12) 

The respective Boltzmann weights p and I — p are given by 

e-/3ff(^i) 1 

where E — U {X2) — U (xi) is the small energy difference between the two minima. Since we are interested in a window 
of size \E\ < T, as T ^ 0, the probability to have two almost degenerate minima situated at xi and X2 with an energy 
difference E can be considered to be uniform in the energy leading to the approximation ^ P{E = 0,xi, X2)dE. Here 
we choose the convention that y = X2 — xi is of arbitrary sign. Then to avoid double counting, we define xi as the 
true minimum, so that E >Q. As in the study of the rare events in j^S, 34], the average over E of the quantities of 
interest can be performed by integrating over E g [0, +oo[. More precisely, we introduce the symbol << .. >> to 
denote the following average over configurations with two degenerate minima 

<< v4(p, 2:1, 2:2) >>= j dxi J dx2 J dE P{E = 0,xi,X2) A ^p = ^ ^^^_p^ ,xi,X2^ (14) 

At this point, it is important to stress here that the measure << . >> is not normalizable, in the sense that 
<< 1 >> diverges, but has only a meaning when the observable A vanishes when only one minimum exists, that is 
for p ^ 1 or for y — > 0. As first examples of the measure defined above, let us compute the lower cumulants of the 
position. The variable {x— < a; >) is {l—p){xi — X2) with probability p and p(a;2 — a;i) with probability 1— _p, leading 
to 



C2(T) = < {x- < X >f > ^« pil - p){xi - X2f » +0{T^) (15) 
= dxi J dx2P{E ^ 0,xi,X2){xi - X2)'^ (16) 



Similarly, the fourth cumulant 



C4(T) = <{x- <x>)^-3<{x- <x >)2 >2> (17) 

can be computed as follows 

C4(r) =« p(l -p){l-6p + 6p^){xi - X2)'^ » +0{T^) = + 0(t2) (18) 

The vanishing result is thus the consequence of a special cancellation of two terms of order T. These cancellations 
actually occur at all orders, as already noted in poj for the random directed polymer model. We now present a proof 
within the effective model (|12|) . 
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B. Generating function of thermal cumulants at order T in temperature 

We now consider the partition function of the effective hamiltonian (|12|l with a source A 

ZeffiX) = J dxe-^'"""'^'' = pe^^^i + (1 - p)e-^^^ (19) 

The generating function of disorder averages of the thermal cumulants can be now expanded in temperature as as 

lnZe//(A) = -XxT+ « In (p + {I - p)e-^(^^-^i)) >> +0{T^) (20) 

Since xi has been defined as the true minimum Xmin of the Hamiltonian H{x), and since the variable y = X2 — xi 
has zero odd moments, we may rewrite 

lnZe//(A) ^-\x—+]^[«\n{p+{\~p)e^y)+\n{p+{\~p)e-^y)»\ (21) 



-Xx:^+- dxi dyP{E = 0,xi,xi+y) dE\n[^ 

^ J ~oo J -oo JO 



(l + e-s/^)2 

It turns out that the integral over the energy difference variable E yields the very simple result 

+°° (l + e-^^'^e^y) (l + e^^/^e'^y) _ dw (l + w'^ + 2w cosh{\y)) _ T ^ ^ 



dElnl^ : -=T Inf^ ^ '-^^ = -X^y^ (22) 

Using the symmetry in y —y, the final result reads 

lnZe//(A) =-Xx^+-X'' dyy^D{y) + OiT^) (23) 

^ Jo 



where 



dxiP{E = 0,xi,xi+y) (24) 

-oo 

represents the probability to have two degenerate minima separated by a distance y > 0. 

The conclusion is thus that the exact identities for the thermal cumulants derived from the statistical 'tilt' 
symmetry [l^ are in fact true at order T for all the models lO that do not present the statistical tilt symmetry. 



C. Disorder averages of the moments at order T in temperature 

Similarly, the even moments of {x— < x >) have the following small-T expansion 

<{x^ <x >)2'= > =« [p{l - pf" + +(1 - p)p^''] {xi - x^f" » +0{T^) (25) 

rp /• + QC 

= -J^ y^^D{y) + 0{T^) (26) 

where D{y) is defined as in H24|l . So the whole distribution D{y) is important to characterize the low temperature 
properties. The comparison with (O shows that there are a lot of cancellations in the disorder averages of the 
cumulants. This phenomenon was already emphasized in j2(J| for the random directed polymer. 



D. 'Chaos' in position at order T in temperature 

Another interesting observables are the 'chaos observables' that contain information about the configurational 
changes induced by temperature shifts. For instance, among the various chaos observables that have been discussed 
for the Sinai model ,25j. the following quantity 



dTi,T2 = {< X >Ti - < X >T2Y 



(27) 
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characterizes the change in the mean thermal position between two temperatures. Again, at order T in temperature, 
this quantity wiU be governed by configurations with two nearly degenerate minima 







rfo,T=«(l-p)na;2-a;i)' »+0(T^) = (21n2-l) / y'D{y)+0{T') (28) 
More generally, higher moments have the following temperature expansion 



4-00 



{<x>T-<x >T=oY^ =« (1 - pf''{^2 - » +0(T2) = Tmk I y^^Diy) + 0{T') (29) 

where 







1 ^2k-l 







(1 + w)^*; 

One can also compute, in the limit where T1/T2 = z is fixed and Ti 



^k = 2 I dwj^-—^ (30) 



r + <yo 

{<x>T,-<x >Tj2fe = TiM,{TjT2) / y^^Diy) (31) 

Jo 



with 



^ dw {w — w^Y^ 

W (1 + W)2fe(l + 



M,(z) = 2/ — ,UfcnL,.m (32) 



In particular, for a small change 5T <^ T, the change in position is characterized by the moments 

i<x>T-<x >T+6T?^ - / dw-——:{\nwf'' / y^'^D{y) (33) 

T 7o (l + w)*'' Jo 

E. Discussion 

All quantities at order T in temperature may thus be exactly computed with the simple eff'ective thermal measure 
((T^ where the statistics of the parameters {xi,X2tP) is defined by ((T^ . For the case of the 'toy model' (quadratic 
well plus Brownian potential), which present the 'statistical tilt symmetry' we have already studied in the 
configurations with two nearly degenerate minima and showed that they give the exact result at order T for the 
generating function of thermal cumulants. We will now consider models which do not present the statistical tilt 
symmetry. 

F. Explicit study of the Brownian potential on finite sample 

In this section, we consider a pure Brownian potential , i.e . Hp = ^ in jSJ, on a finite sample x G [0, i]. Many 
properties of this model have already been studied |23L |26L 1271 123 . l29l| . We first derive the exact thermal cumulants of 
the position, and we then compare with the contribution of configurations presenting two nearly degenerate minima. 

1. Exact temperature series expansion for the thermal cumulants of the position 
The generating function of the disorder averages of thermal cumulants is given by 



ln<e-^^> = InZL(A) - InZL(O) (34) 
where Zl(\) represents the partition function with source A 

Zl[\)= /^dxe-^"e-''^(") (35) 
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The full probability distribution of this partition function is exactly known |28l l29|. For our present purposes, it is 
simpler to work with the following formula in Laplace transform with respect to the length L |29| 



+00 



InZi(A) = In (^^) + ^(1) + ^ [V'^^ + a^-a 



1 + 



2 +v. 1 + 



where ^{x) — with the following notations 



A 

2 fZJ 



(36) 

(37) 
(38) 



we need to perform a series expansion in A, i.e. in a = in formula (|36() . Since the only term odd in A on the right 



For completeness, a derivation of H36(l via path- integrals is given in Appendix El To obtain the thermal cumulants 
hand-side is simply (— a^), we get 

L 



< X > 



(39) 



as expected by symmetry, and all other odd cumulants identically vanish, as a consequence of the symmetry x L — x 
of the random potential V{x). We now consider the even cumulants. The cumulant of order two is given by 



OJ 



+ C30 



dLe 



\< x'^ > - < X >2 



1 



(a/32)2 



(40) 
(41) 



We may now expand in v (|38|l to obtain the small temperature expansion, which also corresponds to the large L 
expansion 



dLe 



'(< x2 > - < X >2) = 



-r{i)+o{v) 



T 



2rp3 



After Laplace inversion, the final result reads 



C2(r) = (< a;2 > - < x >2) = 



T 



47^^5/2 12cr3/2^3/2 o 

/2 _ 7r3/2r3£V2 ^ 



2^5 

2,., 0-5/2^1/2) 



L 



6a3/2 



+o(^^) 



Similarly, the cumulant of order 4 has the following low-temperature expansion 



(42) 
(43) 

(44) 

(45) 
(46) 

(47) 
(48) 

Note that here the next correction is of order T^. We will now show that this result at order T exactly comes from 
the configurations with two degenerate minima. 



C4(r) = (< x'^ > -4 < x^ >< X > -3 < x^ >^ +12 < a;2 >< x -6 < x >"*) 

y3^5/2 T,V2j.5j^3/2 ^7/2^7^1/2 y8 j.9 

-V'""(l)^ + 0(-57^) 



10V7Fcr3/2 12(j5/2 8cr7/2 

More generally, the leading contribution to the cumulant of order {2k) reads : 

(-i)^+ir(fc-i) 



7r(2/c -h l)tT'=-i/2 
The final result is thus that at low temperature 

L TA2 



2-^ ji2fc-l£fc+l/2 



In < e-^^ > = -A- + -^L3/2 _^ o(y3) 
2 6v7r(T 
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2. Study of configurations with two degenerate minima 

It is straightforward to obtain that the probabihty to have two degenerate minima situated at xi < X2 for a 
Brownian potential on [0, L] reads 

P(£;-0,xi,X2) = — -— — ,3,, rj (49) 

The total probability that it happens with a distance y between the two minima reads 

D{y) =e{L-y) [ " dx,P{E = 0,x,,xi+y)^ f (50) 



with the following moments 



dyy-D{y) = —— ^ (51) 



This means that configurations with two degenerate minima appear with a probability of order \j\faL and that the 
distance between the minima is of order L. In particular, the second moment 



^ , Li 
dyy'Diy) = —= (52) 



exactly corresponds to the result H48|l via jnj. 
The relation ()26|l for the moments yields 



<(-^-<^»">-Ij( ^"^(^)-I V^((4fe-l)) 



whereas the relation (|28|l reads 



(2 In 2- 1) 



{<x>T-<x >T=oY = TLi ' ^ ' + O(T^) (54) 



G. Explicit results for the biased Brownian potential on the semi-infinite line 

We now consider (O with Hq = fx with / > on the semi- infinite line x G [0, +00 [ [sj- Again, we compare exact 
results with the contribution of configurations with two nearly degenerate minima. 

1. Exact results for the thermal cumulants 

The distribution of 

ZfiX)= dxe-^"e-'^(-^"+^("» = Zo(A + ^/) (55) 

Jo 



is exactly known [28|, |29|, |30|, |31| , and in particular we have |29|, |3 



lnZ,(A)^ln(^)+-^^-^(l + :^ + ^) (56) 



where ip^x) — r'(a;)/r(x). 

The expansion in A yields the cumulants 



-2 + + — = 72 + -r- 

J O- O" / OCT 

C2(T) ^<a:2>-<x>2 = fif - 1^V/'(1 + — ) = ^ - + 0{T^) (58) 



ci(r) ^<-F> = ^ + — 7/^'(i + ^) = -^ + ^ + o(r3) (57) 

/ cr O" / OCT 



C3(r) = (< a;3 > -3 < a;2 >< a; > +2 < X >3) = — — + — - + O(r^) (59) 

' j4 25s-^ 
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More generally, the leading order at small temperature of the cumulant of order k reads 

1 

Finally, at order T, the exact series expansion of the generating function of the cumulants reads 

a TaX^ 



Ck{T) = -,,^1 + ... (60) 



^p+^ + OiT') (61) 
We will now compare with the analysis of configurations with two nearly degenerate minima. 

2. Configurations with two degenerate minima 
The distribution of the minimum xq for the biased Brownian potential on [0, +oo[ reads 

V 47rCT Jxo X 



In particular, the first moment reads 

p+oo 



dxQXQPy^in{xQ) = — (63) 



in agreement with ci{T) l(5^ at T = 0. 

The probability to have two degenerate minima situated at xi < X2 on [0, L] reads 

P{E = 0, Xi,X2) = Prmn{xi)df{x2 - Xi) (64) 

where the probability df{x) of return to the minimum after a distance x reads 

^ TsbT^'"* '''' 

The total probability that it happens with a distance y between the two minima reads 

Df{y)^ I dxiP{E^O,xux,+y)^df{y)^^=-^e-^ (66) 







In particular, the second moment reads 



dxx^Vf{x) = -J (67) 



2a 



f 



and thus corresponds exactly to the result (|61|l via 
The relation (|26|) for the moments yields 



<(--<- >r > =T ry''D{y) = ^^^^^('4T"'' (68) 



whereas the relation (|28|l reads 

2a 



i<x>T-<x >T=o)2 = T(2 ln2 - 1)-^ + 0{T^) (69) 
H. Conclusion 

We have shown on various explicit examples (the toy model the pure Brownian potential on a finite sample 
in III Fl or the biased Brownian potential in III G|l how the thermal properties equilibrium properties of a particle in 
one-dimensional potentials can be obtained at order T in temperature from the statistical properties of configurations 
with two nearly degenerate ground states. We will now generalize this approach to a model with many degrees of 
freedom. 
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III. LOW-TEMPERATURE PROPERTIES OF THE ID RANDOM-FIELD ISING MODEL 



A. Model and notations 



The Hamiltonian of the one-dimensional random-field Ising chain reads 

N-l N 
i=l i=l 

where the hi are independent random variables with zero mean value hi — and variance 

hl^g (71) 

In the following wc will moreover assume that the distribution of hi is continuous to avoid the presence of an extensive 
number of degenerate ground states. The equilibrium of this model can be studied in details via a Ma-Dasgupta real- 
space RG approach |34] in the universal regime where the Imry-Ma length ,32j . which represents the typical size of 
domains in the ground-state of H70|l 

LiM = » 1 (72) 

9 

is large. A picture of the ground state is given on Figure Q ■ The essential result that will be useful in the following is 
that the join distribution of the length I and the energy gain F = 2\ X^jii \ of Imry-Ma domain reads in Laplace 
transform [84] 

dle-P^PrAF,l) = 0{F > T j)Urj{p)e-^^-^'''>''^-'^P^ (73) 



where 



Urip) 



V 2g 



wr(p) = cothrjf (74) 
V 2g 



The average length of domain walls is then exactly Lim 



1= I dF I dllPrj{F, l)^^=LiM (75) 



i.e. the density of domain walls is simply njjw = ^/Lim- 




FIG. 2: Picture of the ground-state : the zig-zag line represents the renormalization of the random potential V{x) = 2 X^iLi 
seen by the domain walls, where only extrema separated by Al^ > 4J have been kept. Then all maxima are occupied by a 
domain wall B{—\+), and all minima are occupied by a domain wall A{+\ — ). So descending bonds are domains of + spins, 
whereas ascending bonds are domains of — spins 
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B. Statistics of nearly degenerate excitations 



The thermal excitations of nearly vanishing energy are rare events of two kinds : they have been called rare events of 
type (a) and of type (c) respectively in the previous study on aging dynamical properties 01 . These two possibilities 
are the following : 

(i) The excitations of type (a) involve a single domain wall which has two almost degenerate optimal positions : see 
Figure 131 From 113,01 j '^^^ obtain the probability density PaiE = 0,1) of these excitations in the scaling region 
? ~ L/A/ 1 in terms of (|7!^ 



1 /-r /•+00 /■+00 

plarge(^^^Qj-^^__ I ^y^j (Fo , ^1 ) / dhProiTo, 12)5(1 - (h + h)) (76) 



where 7 is a small cut-off. Indeed, as in the discussion of the Brownian potential, the probability of degenerate minima 
diverges at small distance in the continuum (|50|l . but the moments that appear in physical observables will be well 
defined. Indeed, in Laplace transform, we have using H74|l 

/ d/e-f' I p'r^iE = 0, = -—{-dp) / dV^U^^ip) = -—dpurip) (77) 

Jo i^IM Jo ^IM 

(78) 



1 



L 



IM^ J 
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These excitations have thus indeed a length of order I ^ Ljm, but concerns only a small fraction of order l/Fj of the 
domain walls. After Laplace inversion, the length distribution is given by an infinite sum of exponentials 

pr^^(i? = 0, l)dl = J2 n'e-^''''^ (79) 

Here, since we have a discrete model, the probability of degeneracy is regularized by the lattice at small distances. 
Indeed, nearly degenerate excitations of small length I ^ 1 also exist. They are not described by the universal scaling 
regime (|79|l . but depend on specific properties of the initial distribution P(hi). For instance, the excitations of length 
I = 1 correspond to the domain walls that have a neighbor with a small random field hi < T 0. The density of such 
excitations is simply proportional to the density n = I/Ljm of domain-walls in the ground-state, and to the weight 
P{hi = 0) of the initial distribution at the origin 

p,{E = 0,l = l) = -^2P{h, = 0) (80) 

More generally, the statistics of small excitations / = 2,3, .. is governed by the probabilities of returns to the origin 
for a constrained sum of I of random variables. For instance, the density of two-spin excitations reads, assuming the 




FIG. 3: Representation of a two-level excitation of type (a) : a domain wall A of the ground state may have two nearly 
degenerate optimal positions, separated by a distance I if AV = 2^^^^ hi ~ 0. 



B 




A 



FIG. 4: Representation of a two-level excitation of type (c) : a pair (A,B) of neighboring domain walls separated by a distance 
I may appear or annihilate with almost no energy cost if AV = 2^'^-^ ~ 4 J. 



symmetry P{—h) = P{h) 



Pa{E^Q,l = 2) = 



L 



lAI 



dhP^ih) 



(81) 



(ii) The excitations of type (c) involve a pair of domain walls which can appear or annihilate with almost no energy 
cost : see Figure 0] These excitation have by definition a large length I ~ ^ 1 . From |^ , their density 
reads in terms of (fz!^ 



i.e. in Laplace transform, we have more explicitly using (|74|1 



+ 00 



dle^P^PciE = 0,0 



LiM 



IM 



1 - \limP+ ^L]j^jp^ 



(82) 



(83) 



As in H78|) . these excitations of length I ^ Ljm concern only a small fraction of order l/Fj of the domain walls. 
Again, after Laplace inversion, the length distribution is given by an infinite sum of exponentials 



(84) 



In conclusion, there exists on one hand large excitations of length I ^ Ljm S> 1, whose universal probability 
densities are given by the explicit expressions (|79I84|I . and there are on the other hand small excitations of length 
I ~ 1,2,.., whose statistics depend on the initial random field distribution, with the weights l|80l I81|l for Z = 1 and 
I = 2 for instance. 



C. Effective Hamiltonian at order T in temperature 

We note Co — {S'l"''} the true ground state, and we index by a = 1,2,..,Z the possible strings where the spins 
{Si^, ...,Sj^) can be coherently returned with a small amount of energy T 0, as a consequence of the rare 

excitations (a) or (c) described above. These excitations appear with a small probability at low temperature, i.e. 
their are dilute, and the effective hamiltonian at low-temperature is thus given by 

I n'5(^.,^r) I E^^^ifl -5(5., ^fVe-^ fl SiS.,-S^^)] (85) 
ye-R / Q=i ( 1 + e t \ i^i^ J 

where Eq = 7i(Co) is the energy of the ground state configuration and where R is the set of all the spins j that are 
outside all the strings a. 
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D. Specific heat at low temperature 



At low temperature, the effective partition function ()85|l is thus given by a product over independent two-level 
excitations with respect to the ground state 

a 

The free-energy per site reads in the thermodynamic hmit thus reads 



/(T) = -T Jhn^ lln Z^ff = - T £ dEp{E) In (l 



(87) 



where the ground state energy eo = —J — g/{2J) was already discussed in [Sj, and where p{E) represents the 
probability density of excitations of all types and of all sizes 



p(i?) = ^[p„(ii;,0 + Pc(£;,0] 



I 

The linear term in temperature of the specific heat is directly related to the density of excitations at zero-energy 

C{T) = T—p{E = 0) 0(t2) (89) 

E. Edwards-Anderson order parameter 

We now consider the generating function of the thermal fluctuations of a given spin Sk 



ZeffiX) = In < e-^A-s, > (90) 
-™ax (1 -I- e~"^e^-^*^^'" J 

lnZe//(Afe) J2 ^(*" ^ <ia)ln^^ ^ » (91) 

1 + e T 



where << .. >> denotes the average with respect to the properties of the nearly-degenerate strings a. Using now 
5^.°^ = ±1 and the results of ^ for the average over the energies Ea, we get 

r+oo 

In < 6--^'=^'' >^TXl / dl lp{E = 0, + 0(T2) (92) 

The thermal fluctuations of a spin are thus directly proportional to the probability /q^°° dl lp{E — 0, l) that the spin 
belongs to a nearly degenerate excitation. In particular, since S*^ = 1, the Edwards- Anderson order parameter reads 

qEA = < Sk>^ ^l-2T dl lp{E = 0J) + O{T^) (93) 

The vanishing of higher cumulants in (|92|1 yields further identities, the first ones being for instance 
— J 4< ^fe >2 - 1 , , 8^ f+°° 



3 



-0{T'^) = l--T dllp{E = 0,l) + O{T^) (94) 
3 Jo 



30< Sk >4 - 17< Sk>^ + 2 , ^,^2, _ ^ 46^ /■+°° ^, , _ ^ , . 



< Sk >^ = — + 0{T') = 1 T / dl lp{E = 0J) + 0{T') (95) 

15 15 
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F. Susceptibility 

We now consider the global magnetization 

^E^^ (96) 



N 

i 

m 



N 

4=1 



and the generating function of its thermal fluctuations 



Ze//(A) = In < e-^'» > (97) 
The effective Hamiltonian at order T in temperature H85|) yields the following partition function with source A 

■JV 1 -L T C 



Z.,,(A)= e-'^''^— =e-iE-^.^"' n 

Si = ±l a=l 



(^1 + e^-T") 

where = sg7i(S'|°'''") = ±1 with probabilities (±1). We thus obtain 



lnZe//(A)=«^ln^^ _^» (99) 

1 + e T 



where << .. >> denotes the average with respect to the properties of the nearly-degenerate strings a. Using now the 
results of for the average over the energies Ea-, we finally get 



In particular, the average susceptibility reads 



dl fp{E ^0,l) + O{T'^) (100) 



N 



X= ^ m2 > - < m >2) = 2^ dl I'^piE = 0,l) + O{T) (101) 

so the excitations appear here with a weight P. 

The susceptibility (|101|l we have just obtained from the analysis of the thermal fluctuations of the magnetization at 
linear order in T may also be recovered from the analysis of the response to an external field H at zero temperature. 

Indeed, suppose we start from the ground state at iJ = 0, and we turn on a small external field H > 0. The changes 
in the ground state will be localized on the nearly degenerate excitations of small energy E > and length I such that 
the new energy in field allows to decrease the energy Eh = E — 2HI < : the flip of this excitation yields a change 
of m = Hn the magnetization. As a consequence, the total magnetization per spin is given by 

M{H,T = Q)= / dl dEp{E,l)l ~ 2iJ / dlPp{E^O,l) (102) 
Jo Jo Jo 

in agreement with H101|l . as it should to recover the Fluctuation-Dissipation theorem. 



G. Discussion : the role of small/large excitations on various observables 

Since there are two kinds of excitations, the small ones I 1 and the large ones I ~ Ljm, and since various 
observables are governed by the nearly degenerate excitations with various weights involving their lengths, it is now 
interesting to discuss which observables are dominated by small excitations, and which observables are governed by 
large excitations. 

To simplify the discussion, and to compare with the existing rigorous results 2^ . we now consider that the initial 
distribution P{hi) of the random flelds has the following scaling form 




(103) 



15 



where Vix) is a continuous function, such as for instance the exponential distribution = e ^/2 considered in 

[2^ among other cases. The Imry-Ma length then reads 

i/A/ = — = ,+00^ , - (cte)Tp (104) 

9 Jo dxx^V{x) 

The density of large scale excitations scales as (|78l and l55|) 

plarge^j^ ^ 0) = p'^^^'^ {E = 0) + = 0) ^ j-^^ ^ _|_ ^ |! (105) 

On the other hand, the weight of the smallest excitations I = 1 and I = 2 reads H8UI81|) 

p,{E^ 0,1^1) ^J_^EM. (106) 

1 2C°"dxV\x) , , 

Pa{E^Q,l^2) =-^-L_ (107) 

So even if the initial distribution has a hole at zero field 7^(0) = 0, since the prefactor /q^°° dxV^{x) of excitations of 
length 1 = 2 cannot vanish, the density of small excitations has always for scaling 

pr""(i? = o) = -^ — ^ (108) 

In the regime we consider, where the Imry-Ma length is very large L/m ^ 1, the density of small excitations H108f) 
is thus much bigger than the density of large excitations (|105|l 

^.ma/;^^^Q)^^;arse(^^Q) (109) 

As a consequence, the specific heat is governed by the density of these small excitations p08|l 

C{T) = T^piE = 0) + 0{T^) - Tb^ + ... (110) 
6 

where the numerical coefficient b depend upon the form of the initial distribution via the prefactors of the densities 
of small excitations (|107|l . The leading behavior (|110|l is indeed in agreement with the exact results obtained via the 
Dyson-Schmidt method for various models of the form l|103l) . 

Let us now consider the Edwards- Anderson order parameter H93|l : the density of large excitations weighted by 
their lengths has for expression (|78I83() 

dll[p[-^-s^{E = 0,l)+p,{E = 0,l)]= — i^- + -j= — = — (111) 



whereas the contribution of small excitations (|80l I81|l has the same scaling as the density of small excitations (I108fl 

H 

J2 



^ lp,{E = 0,l) = paiE = 0,l = l) + 2pa{E = 0J = 2) + ^ (112) 



Z=l,2. 



In the regime we consider, where the Imry-Ma length is very large 3> 1, the Edwards- Anderson order parameter 
(|93|l is thus dominated by the large excitations contribution (|lll|l and reads 

r + OO rp 

l-qEA^2T dllp{E = 0J) + O{T^)^ — + .. (113) 
Jo ^'J 

Let us now consider the susceptibility HlOlf) : the density of large excitations weighted by the square of their lengths 
has for expression H78I83|I 
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whereas the contribution of small excitations (|80l I81|l has the same scaling as the density of small excitations ()108|l 

l^Pa{E^0,l)^paiE^0,l^l) + Apa{E = 0J = 2) + ---^^ (115) 



1=1.2, 



In the regime we consider, where the Imry-Ma length is very large Ljm S> 1, the disorder averaged susceptibility 
(|101|l is thus dominated by the large excitations contribution (|114|l and reads 



3.9 



This result is in agreement with the following exact result 



l + O 



(116) 



(117) 



obtained via the Dyson-Schmidt method '23| for the case H1U3|) with the symmetric exponential distribution V{x) = 

As a final remark, we may mention the case of discrete distributions for the initial disorder, such as the binary 
case for instance hi = ±h. The physics is completely different, because the nearly degenerate excitations that exist 
for continuous distributions become now exact degeneracies : the entropy thus become finite at zero temperature |^ , 
instead of a unique ground state, but now there is a gap in low-energy excitations |^. 



IV. CONCLUSION 



In this paper, we have shown explicitly on two kinds of ID disordered models how the analysis of the low-temperature 
properties in terms of two-level excitations allows to fully recover the exact results obtained independently. In 
particular, we have obtained that the thermal fluctuations of a particle in ID random potentials are entirely due to 
the rare samples presenting two nearly degenerate minima, whereas the properties of the random field Ising chain 
are controlled by the rare regions that can coherently flip with a vanishing energy cost. We have also explained how 
the identities on thermal cumulants that are valid for systems presenting a statistical tilt symmetry |l7j actually 
hold more generally for various observables at order T in temperature, as a result of the statistical properties of 
two-level excitations. Finally, for the random-fleld Ising chain, our analysis shed light on the influence of small/large 
excitations on various observables : the speciflc heat is dominated by small non-universal excitations that are more 
numerous, whereas the the susceptibility and the Edwards- Anderson order parameter are dominated by universal 
large excitations, whose length is of order Imry-Ma length : these excitations are rarer than the small ones, but 
involve a larger number of spins. 

In conclusion, since the phenomenology of two-level excitations is very often advocated in the field of disordered 
systems, we hope that the present detailed study on ID models is useful to strengthen these ideas. Moreover, since 
the identification of the sample-dependent two-level excitations and the characterization of their statistical properties 
leads to a very clear picture of the low temperature physics, we hope that other disordered systems will be analyzed 
in the same way in the future. 
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APPENDIX A: EXACT RESULTS FOR THE BROWNIAN POTENTIAL ON A FINITE INTERVAL 

1. Basic path-integral 

To study the equilibrium of a particle in a Brownian finite sample, we need the following path-integral 

Giu,l\uo) = r'^^^ VU{v)e-^ioH^r-'^io^'^y^-'"'''-' (Al) 

Ju(0)=uo 
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whose Laplace transform with respect to / has been given in |23| in terms of Bessel functions 



G{u,p\uo) 
G{u,p\uo) 



where = | 



if — OO < U < Uq < +CXD 

if — oo < uq < u < +00 



(A2) 
(A3) 



2. Generating function of the thermal cumulants 

The generating function of the distribution of the partition function with source H35|l may be computed in terms of 
the path-integral l|Aip as follows 



„-Ax-/3!7(x) 



/+CXD i^U{L)—u 2 

du / VU{y)e-^^<> "^^o dye 

-00 J(7(0)=0 

-00 J\/(0)^0 



dueWG'(w,L|0) 
so that in Laplace with respect to L, we get using (|A3|I 



r+00 x2 

dLe-'^^e-i^^(^'> = / dve^''G(v,uj+ — -10) 
4 



where 



with the notations 



Js Z JO Z 



1 fl 

I3 \ cr 



and where a and 1/ have defined in l|38|) . 
To obtain the moments 



(■+00 +00 / \n p+oo 

Jo „^o -^0 



(A4) 
(A5) 
(A6) 

(A7) 
(A8) 

(A9) 

(AlO) 
(All) 



(A12) 



we need to expand (jASJin q, i.e. to expand the function (|A9|) at all orders in s. The properties of Bessel functions 
(differential equation and wronskian) leads to the following differential equation for (|A9I 



Inserting the series expansion 



2a - 1 



1 



1 



+ 00 



(A13) 



(A14) 



into the difTerential equation (|A13|I . one finds a2n+i = = 1//^^ and the recurrence 

(2n — a)2 — 
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oLin = t:^ J- (Aloj 



leading to 



_ r(i-f-f)r(i + f-f) 

"^'^ - MM"r(n+l-f-f)r(n + l + f-f) ^^''^ 



So we obtain the Laplace transform of all moments 

n+l 



(A17) 



in agreement with '2^ and analytic continuation in n — > yields the result (|36|l given in the text. 
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